Bogomol'nyi bound and screw dislocations in a mesoscopic smectic-A 
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The de Gennes free energy functional of an infinite smectic-A liquid crystal at the dual point is 
shown to be topological and to depend only on the number of screw dislocations and the anisotropy. 
This result generalizes the existence of a Bogomol'nyi bound to an anisotropic system. The role 
of the boundary of a finite mesoscopic smectic is to provide a mechanism for the existence of 
thermodynamically stable screw dislocations. We obtain a closed expression for the corresponding 
free energy and a relation between the applied twist and the number of screw dislocations. 

PACS numbers: 61.30. Jf,61.72.Mn 



A successful and striking analogy between different 
phases of a liquid crystal and those of a superconduc- 
tor has been proposed by de Gennes back in 1972 0. 
It provides an example of how geometric and topological 
considerations lead to similar behaviours in two different 
physical systems ■ This analogy has been further used 
by Renn and Lubensky to predict a new phase called the 
Twist Grain Boundary (TGB) phase J3L, whose existence 
has been experimentally confirmed El IB IE El Ell • 
This phase is the liquid crystal analogue of the Abrikosov 
vortex lattice. It consists of screw dislocations that are 
topological defects. The TGB phase lies in between 
the smectic-A (SmA) and the cholesteric phase. In this 
letter, we show that the de Gennes free energy of a 
SmA has an integrability point at which it is topolog- 
ical. This leads to the existence of screw dislocation 
solutions, whose thermodynamic stability is obtained in 
the mesoscopic limit yet to be specified and for the ge- 
ometry of an infinitely long cylinder of radius R. We 
characterize the smectic-A phase by the following set of 
length scales: the transverse and the longitudinal coher- 
ence lengths £j_ and £y, the twist penetration depth At, 
and the wavenumber q s = %L of the density modulation 
of the nematogen molecules along the z-axis of the cylin- 
der. The cholesteric pitch length defined by L p = |p 
is generally larger than Here d is the separation 
between smectic layers and fco is the wavenumber of the 
Frank director in the cholesteric phase. 

We first show that when written in a dimcnsionless 
form, the de Gennes free energy density T of an effective 
two-dimensional system is solely characterized by two di- 
mcnsionless parameters, namely the twist Ginzburg pa- 
rameter n — jj- and the anisotropy parameter a = 

V^k^hQs oc £| | The dual point E3 is defined by 
K = 1/ \/2. At this point and for a large system, the 
de Gennes free energy in the presence of a finite twist is 
a topological quantity that can be written 



T = 2-kI{\ -a 2 ), 



(1) 



mol'nyi bound |l2j for the liquid crystal. However, from 
it appears that there is no mechanism to select the 
number I of defects, and the free energy is minimum for 
I = 0. A way to thermodynamically stabilize solutions 
with I is to consider a mesoscopic finite size smectic, 
just like the unconventional vortex patterns obtained in 
mesoscopic superconductors . A mesoscopic smectic- 
A for the geometry of a cylinder of radius R is charac- 
terized by R <C L p . Moreover, for R > we show 
that for a wide range of parameters ko and a, a meso- 
scopic SmA can accommodate screw dislocations and the 
selection mechanism for the equilibrium value of I =/= is 
provided by the boundary conditions imposed on the to- 
tal twist through the sample, in a way quite similar to 
the Little-Parks setup in superconducting rings. 

We first consider some of the features of the SmA 
phase and the conditions of its description through the 
de Gennes free energy pj. The order parameter that 
characterizes the SmA phase is the first harmonic of the 
density of the nematogen molecules. It can be written 
as ^(r) = |t&|e** where $ = q s n ■ r. The unit vector n 
is the Frank director that specifies the molecular align- 
ment. For a perfect SmA phase we have n = {0,0,1}. 
The smectic planes are defined by z = Id, (I G Z), so 
that the phase $ = 2irl [3, EE Ell- A screw disloca- 
tion is a spiral staircase structure of the smectic planes 
|l4l llfij characterized by a finite winding number. The 
free energy of a SmA is given by the sum of the Landau 
(Fl) and the Frank (Fp) contributions EH which arc, 
respectively: 



Fr, 



Ff 



i-i 2 

(2) 



-K 2 (n.V x n) 2 + K 3 (n x V x n) 2 . 



where the integer / accounts for the number of screw 
dislocations. This expression corresponds to the Bogo- 



The Frank term Fp is the sum of energies associated 
to bend, twist and splay of n. For a perfect SmA, the 
Frank energy vanishes. Hereafter we consider only the 
twist part in Ff. The de Gennes free energy jll llri Il7j 
is obtained as the limit of the total free energy F = 



2 



Fl + Fp , when the Frank director n slightly deviates 
from the z-axis, namely n z w 1 and both n x and n y are 
much smaller than n z . This approximation together with 
the condition that R -C £ p defines the mesoscopic limit. 
The de Gennes free energy is thus given by [l^j 




+ Cj.|(V ± - iq s 5n ± )«£\ 2 + - 2 -(V ± x 5n ± ) 2 , 

where Vj_ = {d x ,d y }. We have introduced the two- 
dimensional vector 5n± = n — z. Since n 2 = 1, at 
the lowest order 8n± lies in the x-y plane. The lengths 
£x and £y over which the SmA order parameter changes 
respectively in the x-y plane and along the z-axis are ob- 
tained by comparing the corresponding gradient terms 
with the first term in the energy functional (3). This 
sets £j_ = yJC±/\r\ and ^ = y/C\\/\r\. The twist pene- 
tration depth is Xt — yj c±*q£\ri Q • ^ n ^ nc mcan n °ld 
approximation, only the coefficient r depends on temper- 
ature 0, 0, 0| , so that the Ginzburg parameter n = j^- 
is temperature independent. By rescaling all the lengths 
in units of s/2\t, namely defining the scaling relations 
^F = F dGl =f\f\ 2 = \V\ 2 , V2X T x = x, V2X T q s = q s 

and Vi = V2ArVx, the de Gennes free energy (3) can 
be rewritten, up to a constant, in the dimensionless form 

T = « a (|/| 2 -l) a + « 3 |/| a + 

+ |(V X - iq s 5n±)f\ 2 + ^(Vi x q s Sn±)'\ (4) 

where both / and 5n± do not depend on z. The 
anisotropy parameter a 2 = 2k 2 q 2 £, 2 = q 2 C\\/C±_ contains 
|r| and it is therefore temperature dependent. The ex- 
pression is analogous to the dimensionless Ginzburg- 
Landau free energy density of a two-dimensional super- 
conductor in a magnetic field [Tsl llflj , except for the ad- 
ditional term a 2 |/| 2 . This term is proportional to 1/d 2 , 
and it accounts for the interaction between the nemato- 
gen density in SmA layers. The equilibrium equations are 
obtained by minimizing the functional Q with respect 
to /* and 5n±. It gives 

\V ± -iq s 6n ± \ 2 f = 2K 2 f(l-\f\ 2 )-a 2 f 
V± x (V ± xq„6n±) = 2j. (5) 

The second relation is the liquid crystal analogue of the 
Maxwell-Ampere equation. It is written in terms of the 
two-dimensional liquid crystal current density, defined 
by j = Jm(/*Vj_/) - \f\ 2 q s Sn±. This current den- 
sity vanishes for a perfect smectic. Using the phase 
of the SmA order parameter, it can also be written as 
3 = |/| 2 (Vj_$ — q s Sn±). In a superconductor, the van- 
ishing of the in-plane circulation of the current density 
around a closed contour T implies the quantization of 



the magnetic flux. Similarly, in a SmA the vanishing of 
the circulation of the two-dimensional current density j 
around a closed contour implies that the smectic planes 
are lifted by an integer multiple of the layer separation 
d. We define the liquid crystal analogue of the London 
fluxoid by 

y Vj$.d?= jf(j^j5+?,*n x ).tff = 27rf, (6) 

where I is an integer. We now make use of the following 
identity which holds in two dimensions only 

\(V ± -iq s Sn±).f\ 2 = \Vf\ 2 +(V ± xj) + \V ± xq s Sn ± \\f\ 2 . 

The scalar operator T> is defined by T> = dx + idy — 
iq s (5n x + iSn y ). We denote by fl the cross-sectional area 
of the cylinder and by dfl its boundary. We set bound- 
ary conditions by imposing that the system behaves as 
a perfect smectic at large distances, namely |/| — > 1 and 
j — > on d£l. Making use of those boundary conditions 
and of the identity JJ, the free energy Q at the dual 
point k = 1/a/2 can be written, up to a constant, as 

/ T = I \ (1 - l/l 2 - « 2 - IVx x q s Sn±\) 2 + W\ 2 
+ {I -a 2 ) <f ( |-7[2 + q s Snx)-dl. (8) 

Jdfl l/l 

The boundary integral results from both the Stokes theo- 
rem and the boundary conditions. Making use of the re- 
lation ©, it is reduced to the constant term 2irl(l — a 2 ). 
Then the free energy © is reduced to the sum of two pos- 
itive terms. It is minimum when each of them vanishes, 
so that the minimum free energy is given by the so-called 
Bogomol'nyi bound £Q) and the corresponding field con- 
figurations JS} arc now solutions of the Bogomol'nyi like 
equations [l2( 

Vf = 0, and |Vx x q s Sn±\ + a 2 = 1 - |/| 2 . (9) 

It is important at this point to notice the role played 
by the nonlinear term in the free energy J3J in order to 
reach the Bogomol'nyi bound. The second equation in 
© expresses the dual relation between the twist and |/|. 
By eliminating the twist between these two equations, 
we obtain for |/| the nonlinear equation V ± Zn|/| 2 = 
2(q' 2 + |/| 2 — 1). For a = 0, it is a Liouville equation which 
is known to admit families of vortex solutions char- 
acterized by the winding number I that appears in the 
free energy JJJ. The existence of a Bogomol'nyi bound 
for the de Gennes free energy of a SmA and the equations 
© are one of the main results of this letter. 

For an infinite system there is no mechanism to select 
the number I of screw dislocations. The boundary of a 
finite mesoscopic SmA introduces such a mechanism that 
determines the value of / corresponding to stable screw 
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dislocation solutions as a function of the twist. In a fi- 
nite system, the order parameter is not equal to 1 at the 
boundary, so that the boundary conditions used to derive 
© are not generally satisfied. Hence, the identification 
of the boundary integral in with the fiuxoid © is no 
longer possible and the free energy cannot be minimized 
just by using the Bogomol'nyi equations JjJJ. To fix the 
boundary conditions for a finite system, we impose that 
the flux of the twist of the SmA through the cylinder 
is equal to the flux of a cholesteric in the same geom- 
etry. To implement it, we consider first the expression 
of the director n in the cholesteric phase 0, 0] given by 
n c h(x) = (0, sinkox, cosk$x) and obtained by minimizing 
the cholesteric free energy F c h = Fp — K 2 kg / (fi • V x n) . 
The flux of the twist through a cylinder of radius R is 
given by 7r R 2 n c h ■ V x n c h = nR 2 k^. And because the 
total flux of the twist of the SmA is 2TrRSn±(R), we have 



5n ± (R) = 6n ± (R)0 



k R 



ttR- 



(10) 



This boundary condition disregards edge dislocations |2ll 
122^ . For the geometry of a cylinder, the current density j 
has only an azimuthal component. We note that on the 
boundary, Vx^lan = l/R, so that for 



■h 



(11) 



the current density je{R) is negative. Moreover, since the 
current density is positive at the core of a screw disloca- 
tion, there exists a circle along which the current density 
vanishes. It allows to split the domain Q into two con- 
centric subdomains il = Qi U Q2, so that the boundary 
dfli is the zero current density line. The system can now 
be divided into a bulk (Qi) and an edge region (fl 2 )- We 
neglect the interaction between the screw dislocation and 
the edge current, so that the total free energy density of 
the system is now given by T = T(VL\) + ^(f^). Since 
the boundary Dili is the zero current line, we can approx- 
imate J 7 (Hi) by the Bogomol'nyi bound Q 0. For the 
edge free energy, we use expression Q integrated over 
f?2- The twist Tir) = Vj_ x q s Sn± enters the system up 
to a distance of order 1 in units of XtV2 Therefore, it 
decreases from the boundary at r = R with a behaviour 
well described by 



T(r) = Vj_ x q s Snj_ 



ch 



R 



(12) 



where is defined by Hll|) . The value at the boundary 
T(R) = (<j> c /j — l)/R results from the expression of the 
total twist 27r<E> c /j as a sum 2-nl + J" n T(r) of the twists 
associated respectively to the dislocations in the subdo- 
main fli and to the contribution of the subdomain il 2 - 
The twist contribution to the free energy in expression 
(jU) is thus given by J Q ±(Vj_xq s 5n.j_) 2 = (l-$ ch ) 2 /4~R. 



Along the same line of arguments, the quantity j/\f\ 2 = 
V ± $-q s Sn ± is such that ((V± X j)/\f\ 2 + T(r)) = 
0. Then, — j/\f\ 2 restricted to the subdomain Q 2 is also 
given by l|12[). At the dual point we have At = £_l/v2, 
so that the amplitude of the nematogen density modu- 
lation / saturates to its equilibrium value (undistorted 
SmA) over a distance \tV2. With the help of these ob- 
servations, we estimate J-(il2) using a variational ansatz, 
namely that the modulus |/| has a constant value /o over 
a ring of width Ata/2, included in £^2, so that the term 
V]_|/||an = [ljj . This leads to an edge energy of the 
form: 



T 



(I - ^chf 



ttR 



AR 
(I 



(13) 



^^/o 2 + (l-/o 2 ) 2 + 2« 2 /o 2 
R 



We then minimize with respect to /o, which yields 2(1 — 
fo) = ^~^''" > — h 2a 2 . The corresponding expression for 
the extremum value of the edge free energy Ijl4(l together 
with the bulk term 0) gives for the total free energy 



FK*^) = 2^(1- a 2 ) + J[|(Z-$ c/l ) 2 + 2i?V] 



R l 2' 



—i(i-*c h y 

4R 



2lfa 2 ] 2 . 



(14) 



This relation constitutes the second main result of this 
letter. In the limit of a radius R 3> At, the quartic term 
becomes negligible compared to the quadratic one, so 
that the ^c^-dependence of the total free energy appears 
as the envelop of a family of parabolae each indexed by 
the integer I. The system chooses its winding number / in 
order to minimize the free energy. This provides a rela- 
tion between the number I of screw dislocations and the 
twist <fr c h given by the integer part I = [$ c h~ a 2 )]. 
The entrance of the (I + l)-th screw dislocation occurs for 
a value of the twist, so that F; = Fj+i. The entrance of 
the first dislocation is thus given by $ c /j = | + |i?(l— a 2 ), 
which according to l|ll|) gives for the corresponding twist 

the expression ko = -^-^ (-£ + ' J ■ The highest 

(c) 

value fcp of the twist that allows the entrance of screw 
dislocations corresponds to the critical value at which the 

(c) I 2 

system becomes cholesteric. It is given by Kq = W 

and it corresponds to $^ = R 2 j\[2. Using this ex- 
pression in (f 1 1 1> - we obtain a bound for the largest num- 
ber of screw dislocations that a mesoscopic sample can 
accomodate. All these results are summarized in Figure 
^ which presents the behaviour of the free energy F; as a 
function of & c h- In deriving the expression l[14|). we have 
implicitely assumed that the parameters oc fco and 
a 2 can be changed independently. Nevertheless, it has 
been shown experimentally [j| that both ko and a 2 de- 
pend on temperature. However, it has been argued that 



4 



1=0 




5 10 



ch 

FIG. 1: Free energy F; of a mesoscopic SmA as a function of 
the twist <& c h as given by \14\/ with R = 8 and a 2 = 0.9. 



of the boundary of a finite size mesoscopic smectic is to 
provide a mechanism for the existence of stable screw 
dislocations. The Bogomol'nyi equations © have been 
obtained without any specific assumption on the nature 
of the twist field. Therefore, they share some similarities 
with those obtained in other situations like mesoscopic 
superconductors fl9j |. Chern-Simons |24| and Yang-Mills 
[25| field theories. These analogies can be pushed further 
towards other anisotropic layered systems like anisotropic 
High-T c sup erconductors or double layer quantum Hall 
systems [26]. 

This research was supported in part by the Israel 
Academy of Sciences and by the Fund for Promotion of 
Research at the Technion. 



kg might also be changed at fixed temperature 0, Hj| . 
Since we do not know the exact temperature dependence 
of fco and a 2 , we have presented in Figure [3 the phase 
diagram of a mesoscopic SmA in the $ c ^-a 2 plane. It 
exhibits finite areas where screw dislocations are ener- 
getically stable. In conclusion, we have shown that the 
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FIG. 2: Phase diagram in the Qch-a 2 plane, of a mesoscopic 
smectic-A with screw dislocations (I = 1, 2, 3) and R = 8. The 
region enclosed by a given contour corresponds to Fi+i < Fi 
and a screw dislocation with winding number I + 1 is ener- 
getically favoured as compared to one with winding number 
I. Thus, in the area (abed), we have a I — 1 dislocation 
while the area (defe) corresponds to I = 2. The critical 

twist $^ = R I\f2 (see text) at which the system becomes 
cholesteric is = 45.26. 

de Gennes free energy of a smectic-A is characterized by 
two dimensionless parameters k and the anisotropy a. 
At the dual point k = 1/V%, the free energy is topolog- 
ical and it saturates the Bogomol'nyi bound. The role 
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